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ABSTRACT 

We describe the X-ray analysis procedure of the on-going Chandra Multiwavelength Plane (ChaMPlane) 
survey and report the initial results from the analysis of 15 selected anti-Galactic center observations (90° < 
/ < 270°). We describe the X-ray analysis procedures for ChaMPlane using custom-developed analysis tools 
appropriate for Galactic sources but also of general use: optimum photometry in crowded fields using advanced 
techniques for overlapping sources, rigorous astrometry and 95% error circles for combining X-ray images or 
matching to optical/IR images, and application of quantile analysis for spectral analysis of faint sources. We 
apply these techniques to 15 anti-Galactic center observations (of 14 distinct fields) in which we have detected 
921 X-ray point sources. We present logA^-logS' distributions and quantile analysis to show that in the hard band 
(2-8 keV) active galactic nuclei dominate the sources. Complete analysis of all ChaMPlane anti-Galactic 
center fields will be given in a subsequent paper, followed by papers on sources in the Galactic center and 
Bulge regions. 

Subject headings: X-ray source. Galactic survey 
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1. EVITRODUCTION 

The on-going Chandra Multiwavelength Plane (ChaM- 
Plane) survey is designed to constrain the low-luminosity 
accretion source content an d stellar-coronal sourc e luminos- 
ity functions in the Galaxy ( iGrindlav et al. 2003, 2005). To 
achieve this goal we search for low-luminosity Galactic X- 
ray sources in Chandra archival data and attempt to identify 
them using follow-up optical and/or infrared (IR) imaging and 
spectroscopic observations. 

This paper describes the X-ray analysis of the ChaMPlane 
survey. To illustrate the procedures, we describe the initial X- 
ray results from the analysis of 15 selected anti-Galactic cen- 
ter (anti-GC) observations (90° < I < 270°; Table[lli. |2lde- 
scribes the custom analysis tools developed for the project. ^ 
presents initial results including logA^-logS' distributions of the 
sources and basic source properties using a new X-ray spec- 
tral c lassification technique - quantile analysis jHong et alJ 
l2()0l . 

In subsequent papers, we will present the X-ray results of 
all anti-GC fields and the Galactic-Bulge and Center region 
sources. The overview of the project is given in Grindlav et alJ 
(|2005). The methodology an d initial resu l ts of the optical and 
IR surveys are presented in "Zhao et alj J2005h . iRoeel etaH 
ll2005 ) and Lavcock et al. ( 2005) , respectively. 



2. DATA ANALYSIS 

The overall data analysis procedures can be grouped into 
two stages. First, we search for X-ray point sources using 
a wavelet detection al gori thm (wavdetect; Mallat. S. 1998; 
iFreeman et alJl2002l ^2.11 . Second, we determine various 
source properties using simple aperture photometry ( ^2.2> . To 
apply the above algorithms consistently on the large ChaM- 
Plane data set, we employed a custom X-ray analysis tool 
(XPIPE) and have developed post-XPIPE procedures (PXP). 
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Both tools are primarily based on CIAO tools (version 3.1).^ 
XPIPE is an X-ray analysis tool developed for the Chandra 
Multiwavelength Project (ChaMP) which is a high-latitude 
extra-Galactic survey for active galactic nuclei (AGN) (Kim 
et al. 2004a and 2004b; K04 hereafter). The scientific goals 
of the two projects, ChaMP and ChaMPlane, are as different 
as the survey regions - Galactic Plane survey versus high- 
latitude survey. However, both projects share similar require- 
ments for the analysis, namely, searching for faint X-ray point 
sources in Chandra ACIS archival data. This analysis similar- 
ity and the success of XPIPE in the ChaMP analysis (K04) led 
us to adopt XPIPE as our primary tool. The detailed descrip- 
tion of XPIPE can be found in K04, and in what follows we 
only briefly review the basics that are relevant for ChaMPlane. 

In addition to XPIPE and in order to meet ChaMPlane- 
specific analysis requirements, we have also developed our 
own analysis tool, PXP, which is described below. PXP uses 
XPIPE outputs and performs analysis that is optimized for 
Galactic sources but also of general use. 

2.1. Source detection by wavdetect 

For each observation, we start the analysis by using XPIPE 
on the level 2 data products generated by the Chandra X- 
ray Center (CXC) standard data processing. First, XPIPE 
removes the residual artifacts that may be present after the 
CXC standard data processing. XPIPE then selects events in 
the good time intervals (GTIs)-^ which are selected to have 
fluctuations < 3(t from the mean background rate. 

For source detection, XPIPE uses wavdetect with a signif- 
icance threshold of lO"*" (about 1 possibly spurious detection 
per CCD, see K04 for the detailed analysis of source detec- 
tion efficiency with wavdetect in XPIPE) and a scale parame- 
ter varying in seven steps between 1 and 64 pixels to cover a 
wide range of source sizes. XPIPE applies wavdetect in three 
separate energy bands, Sx, Hx and Bx (Table |2li, using the 
exposure map generated at 1.5 keV. Note that the wavdetect 

^ http://cxc.harvai'd.edu/ciao. 

^ These are different from the GTIs set by the CXC processing (K04). 
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TABLE 1 

15 SELECTED CHAMPLANE OBSERVATIONS IN THE ANTI-GC REGION (SORTED BY GALACTIC LONGITUDE EXCEPT FOR THE STACKED DATA) 



Obs. ID 


Target 


/ 

(deg) 


b 

(deg) 


Exposure" 
(ksec) 


(xlO-^cm^^) 


No. of 
Sources'^ 


Obs. Date 
(y-m-d) 


CCDs Aim'' 
used ACIS- 


2787 


PSRJ2229+6114 


106.64900 


2.94848 


91.5 


1.04/0.40 


91 


02-03-15 


0123.. 67.. I 


1948 


3EG J2227+6122 


106.64901 


2.94965 


14.7 


1.02/0.40 


20 


01-02-14 


0123.. 67.. I 


755 


B2224+65 


108.63800 


6.84522 


47.6 


0.44 / 0.23 


78 


00-10-21 


..23.5678. S 


2810 


G116.9+0.2 


116.94330 


0.18420 


48.8 


0.48 / 0.28 


97 


02-09-14 


0123.. 67.. I 


2802 


G127.1+0.5 


127.11343 


0.53889 


19.2 


0.89/0.31 


44 


02-09-14 


0123.. 67.. I 



782 


NGC1569 


143.68323 


11.24151 


93.3 


0.42/0.11 


89 


00-04-11 


0123 


.5.7.. S 


650 


GK Persei 


150.95713 


-10.10413 


90.3 


0.20/0.10 


83 


00-02-10 


. .23 


.5678. S 


2218 


3C 129 


160.42891 


0.13717 


30.2 


0.63 / 0.24 


25 


00-12-09 


....456789 S 


676 


GRO J0422-^32 


165.88229 


-11.91286 


18.7 


0.19/0.09 


67 


00-12-09 


0123 


..6.8. I 


2803 


G166.0-^4.2 


166.13140 


4.33778 


29.3 


0.38/0.19 


53 


02-01-30 


0123 


..67.. I 


829 


3C 123 


170.58315 


-11.66140 


46.3 


0.59/0.12 


52 


00-03-21 


. .23 


.5678. S 


2796 


PSR J0538-^2817 


179.71974 


-1.68585 


19.4 


0.80 / 0.26 


31 


02-02-07 


. .23 


.5678. S 


95 


A0620-00 


209.95777 


-6.54014 


41.2 


0.29/0.18 


61 


00-02-29 


. . .3 


.5678. S 


2553 


Maddalena's Cloud 


216.73098 


-2.60034 


24.5 


0.93/0.27 


59 


02-02-08 


0123 


..67.. I 


2545 


Ml-16 


226.80033 


5.62592 


48.6 


0.13/0.12 


71 


02-02-11 


. .23 


.5678. S 


=52787 


2787 & 1948 Stacked 


106.64900 


2.94848 


106.2 


1.04/0.40 


90 




0123 


I 



''Sum of the goo d time intervals (GTIs . |2Ji; 

''Two estimates; Schleael et al.' ("1998^ / iDrimmel et alJ ^2003^ . The total sum along the line of the sight, averaged over source positions. 
§2.2.2 for the differences in the two models. 

^Number of valid (level 1) sources found in the Bx band. See Tables 2 and 4 for the definition of bands and levels. 
''Aim point detector: I = CCD 3 and S = CCD 7. 

''The combined data (CCD 0, 1, 2, 3) of two stackable observations (1948 and 2787). See jOand Table|3l 



See 



routine in CIAO version 3.1 includes source position refine- 
ment procedures that had to be applied separately in the old 
version of XPIPE based on CIAO version 2.3 (K04). 

When multiple observations of a similar region of the sky 
are available, it is possible to combine the data sets to look for 
fainter sources and to study variability. Particularly near or 
at the GC region, ma ny stackable observations are available 
JGrindlav et al J 120051) . We have devised PXP so that it can 
employ wavdetect on stacked images. We consider multiple 
observations stackable if their aimpoints are within 1' of each 
other and the aimpoint detectors are the same type (ACIS-I 
or ACIS-S). Among the stackable observations, we designate 
the one with the longest exposure to be the base observation 
of the stacked image. For the analysis convenience, we as- 
sign a separate Obs. ID for the stacked data {xyyyy, where x 
is > 5 and yyyy is the base Obs. ID) and we assign the CCD 



TABLE 2 

Energy bands for the ChaMPlane X-ray analysis 



Type 




Name 




Range 


XPIPE wavdetect band 












Soft 


Sx 


0.3 


- 2.5 keV 




Hard 


Hx 


2.5 


- 8.0 keV 




Broad 


Bx 


0.3 


- 8.0 keV 


Conventional band 












Soft 


Sc 


0.5 


- 2.0 keV 




Hard 


He 


2.0 


- 8.0 keV 




Broad 


Be 


0.5 


- 8.0 keV 



ID of the stacked data to be 3 for ACIS-I and 7 for ACIS- 
S observations. For stacking, we only use the data of CCD 
0, 1, 2, 3 for ACIS-I and the data of CCD 7 for ACIS-S ob- 
servations. For the aimpoint, detector response function, and 
other required parameters of the stacked data, we employ the 
same of the base observation. In the anti-GC fields (Tabled, 
Obs. ID 52787 is the stacked data of two observations (1948 
and 2787) and Obs. ID 2787 is the base observation of the 
two. 

Before we st ack the data, we apply astrometric corrections 
(see also ^2.3> . First, we calculate the coiTection for the as- 
pect offset of each observation."^ Second, we calculate the 
boresight offset of each observation relative to the base obser- 
vation using the aspect-coiTected position of the wavdetected 
sources. The final boresight offset of each observation is de- 
rived by an iterative procedure using a relatively more reli- 
able subset of matching source pairs. The iterative technique 
is identical to the boresight correction procedure employed 
betwee n X-ray and op tical data for the optical ChaMPlane 
survey dZhao et alJ2005D except that the procedure is applied 
between X-ray and X-r ay data using the estimated 95% X-ray 
positional eiTors ( ^2.3t . 

For the purpose of stacking the X-ray data pointed to a sim- 
ilar region of the sky (i.e. large overlap in the field of view), 
direct comparison of the boresight of an X-ray observation 
to another X-ray observation is more efficient than compar- 
ing it to observations (or catalogs) in another wavelength be- 
cause the number of matching pairs found in the former case 
is likely larger than that in the latter and a large number of 
matching pairs usually provide a reliable estimation of the 

http://cxc.harvai'd.edu/ciao/threads/arcsec_correction 
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TABLE 3 

Source detection in the stacked image (Obs. ID 52787) of two 

OBSERVATIONS (1948 AND 2787) 



Obs. 


Exp. 


Offset (R.A., Dec.) 


"No. of 


Common 


ID 


(ksec) 


Aspect Boresight'^ 


sources 


1948 2787 


1948 


14.7 


-1.82", 0.62" -0.43", -0.05" 


20 


- ''18 


2787 


91.5 


-0.40", 0.03" 


85 


■'18 


52787 


106.2 




90 


20 80 



''Number of valid (level 1) sources in CCD 0, 1, 2, and 3. Note that 
Obs. 2787 has 6 sources in CCD 6 and 7 (see Table0. 
''Number of common sources found in bo th data set. This is deter- 
mined by the X-ray positional error ( 82. 3t . 

^Relative to the base Obs. ID 2787 after correcting aspect offsets. 
Among 18 matching pairs, 12 pairs are selected for calculating the 
final boresight offset. 

''The other two sources of Obs. ID 1948 are located in the chip gap 
of Obs. ID 2787. 



relative boresight offset. In the case of Obs. ID 1948 and 
2787, we have 18 matching pairs out of possible 20. Accord- 
ing to so urce n et count (> 10), off-axis angle (< 8') and level 
(= 1 ; see ^3.1> . 12 of them are selected for calculating the final 
boresigh t off set at the end of the iteration procedure (Table |3] 
and see ^2.3i . 

Once the aspect and boresight offsets are calculated, PXP 
reprojects all the stackable data (XPIPE-screened event files) 
onto a common projection point (the aimpoint of the base ob- 
servation) accordingly, and applies wavdetect on the stacked 
image in the Bx band using the total exposure map (at 1.5 
keV) for the stacked data. In the case of Obs. ID 52787 with 
respect to two Obs. IDs 1948 and 2787, we detect 8 new 
sources and miss 5 sources by stacking (TableO. 

2.2. Source properties by aperture photometry 

Subsequently, PXP employs the source detection results of 
XPIPE (or PXP for stacked data) in the Bx band and extracts 
basic source properties via aperture photometry on XPIPE- 
screened event files using the "conventional" bands - Sc, He 
and Be- shown in Table |2] Note that while XPIPE does de- 
tection in three bands, we only consider detections in the Bx 
band. XPIPE also performs an aperture photometry using all 
6 bands in Table|2] but we implement a separate aperture pho- 
tometry in PXP (modelled after the one in XPIPE) to be effec- 
tive in crowded fields (e.g. GC) and to provide more versatile 
outputs that are particularly useful for describing the diver- 
sity of source populations found in the Galactic Plane. The 
aperture photometry in PXP is also designed to be compatible 
with the stacked data. 

Note also that wavdetect provides basic source properties 
(e.g., net counts). However, we employ wavdetect with a 
39% inclusion radius of the point spread function (PSF; r39%: 
39% of source photo ns lie within the circle^), which is recom- 
manded by Freeman et alJ J2002t for the simple algorithm of 
source characterization in wavdetect. Therefore, source prop- 
erties reported by wavdetect may not be accurate due to the 
large missing fraction of the source region particularly when 
the source is near large regions of diffuse emission. 

' The rj-% values are calculated at 1.5 keV from the Chandra calibration 
data psfsize_20010416.fits i http://cxc.harvard.edu/cal/Acis/ 1. 



Table0]summarizes the key parameters of the aperture pho- 
tometry in PXP. We define the basic source region using a 
circle around the source position (< r95%); the background re- 
gion is defined by an annulus (2 r95% < r < 5 r95%). For practi- 
cal purposes, we limit the range of r95% to 2.5" < r95% < 40".*' 
If the background annulus overlaps with neighboring source 
regions we exclude the neighboring source regions from the 
background region (Tablel^. 

Source regions that overlap are handled in the following 
manner We divide the source region into a core (circle) and 
a shell (annulus). We refine the source region to be the sum 
of the core and a pie sector of the shell that excludes the com- 
mon sector with the neighbor's source region. The core size 
is determined empirically to include as many source photons 
as possible, while minimizing contamination from neighbors. 
If the source region overlaps with multiple neighbors, the 
core radius is determined by the nearest neighbor, and the pie 
sector excludes all the common sectors with the neighbors' 
source regions. 

Fig-fflshows examples of overlapping source regions. The 
shaded region indicates the refined source region: the core 
and the uncommon pie sector (yellow shade for one and blue 
shade for the other source). In the left panel of Fig. ^ where 
the overlap is relatively small, the core region does not over- 
lap with the neighbor's source region. In the right panel, the 
core radius is set to be one-third of the distance between two 
sources due to relatively large overlap. Simple aperture pho- 
tometry using the source region with r < rgs% produces al- 
most identical spectral properties for the two sources in the 
right panel; XS00782B7_001 and XS00782B7_002.^ Aper- 
ture photometry employing the refined source region reveals a 
si gnific ant difference in the spectral types of these two sources 
( j2T3l andFig.l4l. 

Note that the above correction for overlapping source re- 
gion is modified from the original correction implemented 
in XPIPE (K04; Kim 2005) in order to al low re latively high 
source counts and easy quantile analysis (' ^2.2.3t . 

XPIPE sets various flags on each source detected by wavde- 




''5<S00782B7_016 



4180 4185 4190 4195 

X (pixel) 




XS00782B7_001 



3700 3705 3710 3715 

X (pixel) 



Fig. 1. — Examples of overlapping source regions; relatively small over- 
lap (left) and large overlap (right). The refined source region (yellow shade 
for one and blue for the other source) consists of a core and the uncommon 
pie sector. In the case of XS00782B7_001 and XS00782B_002 (right panel), 
this refinement in aperture photometry reveals these two sources have sig- 
nificantly different spectral types, which are otherwise indistinguishable (see 
i|2.2.3i . See the electronic ApJ for the color version of the figure 

* The actual minimum value of rg^% is 2". In order to compensate for a 
possible positional error, we add 0.5", which is the expected posit ional eiTor 
for ~ 20 count sources at < 3' away from the aimpoint. See i|2.3l 

^ Source ID follows the ChaMP convemtion: XS{Obs. ID) {energy 
band} {CCD ID }_{ Source No. ) and the Bx band is the B band in ChaMP. 
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TABLE 4 

Aperture photometry in ChaMPlane analysis 



Source 
Region Overlap 



Condition 



Core Radius 

('■.) 



Refined 
Source Region 



Background 
Region 



No 



Yes 



r < Fc 



■ '■95% 



A < l.5r'g„, 



A > ^68% + rgy/^ 



^ ''95% 



A < r6g% + r'g„^ 



A/3 



r < rc and 
pie sector in 

rc <r< rqs% 



2 r95% < r < 5 r95% and 
r" > rg'j^j, for all neighbors 



Notes. — 2.5" < rq5% < 40". A is the distance between the source and the nearest neighbor, and r^j^ is the 95% PSF radius of the nearest 
neighbor. At a given position in the sky, r is the distance from the source and r" is the distance from neighbors with the 95% PSF rg'so/^. Note 
that when there is an overlap, r'gje/^ ~ rg5% because of the relatively small change of the PSF size compared to the change in source position 
(< 2%). When overlapping is severe (rc < ^39% and 1 /R' < 0.3), we return to a simple aperture photemetry (flag=142). See Eq. J2j for R' . The 
PSF radii are calculated at 1.5 keV. 



tect depending on source properties such as proximity to a 
chip boundary. These flags allow us to determine the validity 
of each source. A subset of the flags are set manually by the 
V&V step (K04) in order to remove residual artificial sources 
such as those stemming from frame transfer streaks during 
read-outs. PXP inherits the XPIPE V&V-based source flag 
information and sets additional flags relevant to the aperture 
photometry routine in PXP. 

Table |5] lists the flags and their definitions. For example, 
when the source region overlaps with neighbors' source re- 
gions, we flag the source (141). If the core for the refined 
source region is very small [rc < ^39%) and the area of the 
uncommon pie sector is less than 30% of the shell area, we 
return to the simple aperture photometry because the refined 
aperture photometry is also subject to a large uncertainty and 
we set another flag (142) on the source for further analysis. 

For the initial analysis, we extract three groups of source 
properties from the aperture photometry: source count and 
rates, flux, and quantiles. 

2.2.1. Net counts and net count rate 

For a given band, the number of net counts (A^net) of a source 
without source region overlap is derived from the number of 
counts (A^src) in the source region (with area As = Trr^j^J, the 
number of counts (A^bkg) in the background region (with area 
Ab), and the relative ratio {R) of the two areas scaled by the 
exposure map (e^) values of the regions. 

A^„et = A^src-/?A^bkg, and 7? = ^^,^^^, (1) 

where e^lpeA is the mean value of Cp in A. Note that ideally 
A 'e]j\p^A = J2peA ^P' when A is small (often true for source 
regions), the former produces more reliable estimates of the 
exposure sum because exposure maps are usually pixellated. 

In the case of source region overlap, we assume azimuthal 
symmetry of the PSF and we derive A^sic from the number of 
counts (A^core) in the core and the number of counts (A^pie) in 
the uncommon pie sector (with area Ap), 

A^sic = A^coi-e + R' A^p,e , and R' = ^xl^ ^ (2) 

Ap 

where A'p is the area of the common pie sector and we use 
Eq. Q for A^net- 

The net count rate (rnet) is defined by the ratio of the number 
of net counts to the exposure time (T) and it is scaled by the 



ratio of the mean value of exposure map within the source 
region to the value at the reference point, max(ep), 

_ A^net max(ep) 

'"net -^7^^ ■ {j) 

i ep\p£As 

The max(ep) is the maximum of the exposure map value 
of the chip. This exposure map scaling is implemented for 
easy conversion from count rate to flux regardless of source 
position on the chip. 

PXP calculates the number of net counts, the net count rate 
and their errors for each source in the Sc, He, Be, and Bx 
bands (Table|3- 

2.2.2. Flux 

To estimate the X-ray flux from a source, PXP uses sherpa^ 
to calculate the count rate to flux conversion factor. To do 
so, we need an accurate detector response function, a column 
density (A^h) to the source, and a reliable spectral model for 
X-ray emission. 

The detector response function varies with source position 
on the chip. However, we simply use the redistribution matrix 
file (RMF) and the auxiliary response file (ARF) calculated 
at the reference point (exposure maximum) for each CCD be- 
cause spatial variation of RMFs and ARFs is relatively small 
(usually < 5%).^ Note that the ARFs generated by mkarf ac- 
count for the known temporal degradation of the low energy 
efficiency of ACIS detectors. '° For stacked data, we simply 
use the RMFs and ARFs of the base observation for the initial 
estimation. 

For the spectral model or column density, as we do not 
have an a priori prescription or estimate for most detected 
sources, we use several simple models to estimate the source 
flux, viz. power-law models with Fph =1.7, 1.4 and 1.2 likely 
to be appropriate for accretion powered sources such as back- 
ground AGN and X-ray binaries. We also use Raymond- 
Smith, MeKaL, thermal Bremsstrahlung, and Black Body - 
all with kT = 1.0 keV - for thermal sources (e.g., coronal 
sources)." 

http://cxc. harvard. edu/sherpa/thi'eads/index. html 
' The CIAO version 3.2 can generate the spatially variant RMF and ARFs. 
The analysis presented in this paper is performed under the CIAO version 
3.1. 

Implemented in the mkarf of the CIAO tools (version 3.0 or higher). See 
http://cxc.harvard.edu/cal/Acis/ 

' ' Corresponding model names in sherpa are xspowerlaw, xsraymond, xs- 
mekal, xsbremss and xsbbody with xswabs for the column density. 



5 



TABLE 5 

Source Flags used in ChaMPlane Analysis 



Type 


Flag 


Definition 




111 


False source by hot pixels or by bad bias values 


False sources by 


112 


False source by bad columns 


113 


False source due to readout streaks by a very strong source 


114 


False source by the FEP 0/3 problem'' 




115 


Double source detected by the PSF effect*^ 




121 


Other spurious sources 


Valid sources but 
source properties 
may be subject to a 


141 
142 


Source region overlaps with neighbors' source regions 

No refinement on source region due to severe source region overlap 


143 


Background region overlaps with neighbors' source regions 


d 

large uncertainty 


146= 


Source at the chip boundary 




147 


Pile-up candidate (> 100 cts/ksec) 




148 


Uncertain source position by flag =115*^ 




151 


Source is extended 


Other cases 


153 


Target of the observation 


154 


Sources in the target region (e.g. point sources in a target cluster) 




157' 


The same source is found in other observations 




158' 


The same source is found in the stacked image 



Notes. — XPIPE assigns a very similar set of flags (see Table 3 in K04), but we redefine a full set of flags for PXP because many flags depend 
on the details of the aperture photometry. To avoid confusion with the XPIPE flags, the PXP flag IDs > 100 and the XPIPE flag IDs < 100. 
"Provided by the XPIPE Validate and Verify (V&V) record. 
I^ttpi/Zcxc. harvard. edu/ciao/caveatS| 
'^See Fig. 7 in K04. 

''These flags are set automatically. The 95% PSF size for source region is based on psfsize_20010416.fits ( 92.21 . 
■^When the source region contains any pixel that has 10% or less of the maximum exposure value of the CCD. 
'These flags are set automatically, based on Eq. j6|. 



For the column density A^h. we employ tw o models; 
iSchleael et al. ( 1998) and Drimm el etalT ( 1200311 ). In gen- 
eral, .Schlegel et al, (,1998) provide realistic estimates of to- 
tal Galactic extinction, but their estimation near the Galactic 
Plane {\b\ < 5°) is not as reliable as in high-latitude fields. 
The model may overestimate the extinction due to incom- 
plete correction for p oint source contributions in the Plane. 
iDrimmel et al.l J2003!) provide a 3-D model of extinction in 
the Galaxy, which is useful for Galactic sources with known 
distances.'^ However, their model may underestimate the to- 
tal column density since it does not include contributions from 
all known molecular clouds. 

As shown in Table 1 , the A^h estimates by ISchlegel et al.l 
([^98) for the selected anti-GC fields are indeed larger (by 
^ 10-200% depending on the fields) than the estimates by 
iDrimmel et al. (2003). According to Drimmel et al. (2003), 
their model is normalized so that the two models agree well 
for high latitude fields, but as indicated by large disagree- 
ments, both models suffer relatively large uncertainties at low 
latitude fields. 

In addition to intrinsic uncertainties in the models, we do 
not know the distance of most of ChaMPlane sources. There- 
fore, PXP calculates the A^h value at each source position for 
both models, assuming the total column density along the line 
of the sight. The angular resolutions of the models are ^ 6' 
for Schleael et al. ( 1998) and - 20' for Drimmel et al. (2003), 
and A^H is then interpolated on this grid to the source position. 

There is a distance parameter in the model by Drimmel et al. i2003). 
The model has only Galactic components, so that its estimate does not change 
for distances greater than ~ 20 or 30 kpc. 



We ii se Nh/E(B- V) = 5.8 x 10^' cm'^/mag dBohlin et alJ 
[T978I) to convert Sch legel et af ( 1998) values of E(B- V) and 
Na/AiV) = 1.79 X 10^' cm'^/mag (Predehl & Schmitt 199^ 
with the "rescaled" option to account for the small angular 
scale dust structure in using the Drimme l et alj (l2003i) maps. 

For each source, by default, PXP estimates the X-ray flux 
in the Sc, He, Be and Bx band using all seven spectral mod- 
els and two extinction models. PXP is designed to allow easy 
addition of new spectral or extinction models and it can also 
assign a separate spectral and extinction model unique to each 
source. In future, for sources with information on spectral 
type, distance or column density, the X-ray flux will be re- 
vised accordingly. For example, one can estimate the column 
density to a source from the quantile analysis of the X-ray 
spectrum if one can assign a reliable spectral model to the 
source (see ^3.3> . 

2.2.3. Quantiles 

The survey strategy ultimately relies on the optical or 
IR spectroscopy to identify the nature of low luminosity 
X-ray sources discovered from the Chandra archival data 
(Grindlav et al. 2005). However, even at the R ~ 24.5 optical 
ChaMPlane limits (Zhao et al 2005), we do not find counter- 
parts for about half the ChaMPlane X-ray sources (e.g. Ta- 
ble0. In these cases, we rely on the X-ray data for a clue to 
the source nature. Since the majority of the sources possess 
only ^10-20 counts it is important to use a technique that 
is sensitive to low-count statistics and with minimum count- 
dependent bias. 

In the ChaMP project, source classification by the X-ray 



6 



data relies on the X-ray hardness ratio and X-ray colors (K04). 
These quantities are based on a particular choice of three sub- 
energy bands (0.3-0.9, 0.9-2.5, and 2.5-8.0 keV). In ChaM- 
Plane, it became quickly clear that this band choice is not op- 
timal for many highly-absorbed Galactic sources or very soft 
coronal sources due to count-dependent s election effects in - 
trinsic to this particular choice of bands jHong et al J 120041) . 
In fact, for X-ray hardness ratio or colors, there is no single 
set of energy bands that can effectively describe the diverse 
classes of X-ray sources found in Galactic Plane fields. 

Therefore, for ChaMPlane, we employ a new spectral clas- 
sification technique, quantile analysis, to acquire more ver- 
satile measures of X-ray characteristics of sources. First in- 
troduced by Hong et al. (2004), the quantile analysis employs 
various quantiles (median, terciles, quartiles, etc) of spectral 
distributions to reveal and classify various spectral features 
and shapes. Hons et al. (2004) illustrate the technique using 
a quantile-based color-color diagram (QCCD), where the me- 
dian of the distribution is used to indicate the overall hard- 
ness and the quartile ratio is used to classify the general shape 
(concave-up and -down) of the spectrum. 

The quantile analysis does not require sub-division (bin- 
ning) of the energy range, and it takes full advantage of energy 
resolution of the instruments. Therefore, quantile analysis is 
free of any selection effects inherent in the conventional ap- 
proaches, and source classification by QCCDs is uniformly 
more se nsitive th an that by conventional X-ray color-color di- 
agrams dHong et al.^2004.) . 

For a given source, PXP feeds the lists of photon ener- 
gies in the source and background regions together with the 
weighted ratio of the two regions (R) into the routine quan- 
tile.pl (Hong et al. 2004)- '"^ In the case of source region over- 
lap, we use events in the refined source region (the core plus 
the uncommon pie sector) and the area ratio R is given by 



R = 



(4) 



where Ac is the area of the core and Ap the area of the un- 
common pie sector For example, without refining the source 
region, the two sources in the right panel of Fig.^show almost 
identical spectral properties due to source photon mixing; the 
median energy £50% = 3.96(11) keV for XS00782B7_001 
and £50% = 3.96(09) keV for XS00782B7_002. How- 
ever, the above-refined aperture photometry reveals their 
spectral type is indeed quite different; £50% = 4.29(15) 
keV for XS00782B7_001 and £50% = 3.73(12) keV for 
XS00782B7_002 [diamonds in Fig. 111(c)] 

For each source, PXP generates five quantiles, median (m = 
Q50), terciles (233 and Q^j), quartiles (Q25 and Qjs), and their 
errors in two broad energy bands, Bx and Be- Examples of 
qua ntile analysis on the anti-GC fields are provided in section 

2.3. Positional uncertainty and cross-correlating X-ray 
sources 

The wavdetect routine provides the sky position and po- 
sitional errors of each source, but wavdetect underestimates 
the positional errors because the inputs to wavdetect are not 
ideal. For example, we use symmetric PSFs for wavdetect 
but the PSFs become asymmetric at large off-axis angles. In 



Quantile analysis routines are available at 

://hea-www.harvard.edu/ChaMPlane/quantile/ 



addition, the exposure map is calculated at 1.5 keV for prac- 
tical purposes (no information about the source spectrum is 
available in advance; the effective area as a function of en- 
ergy peaks at 1.5 keV) but no ChaMPlane source is expected 
to exhibit purely 1.5 keV emission. For optical or IR identifi- 
cation of X-ray sources, we need a more reliable estimate for 
the uncertainty of wavdetected source positions. To do so, we 
have performed SAGS AC''* & MARX'^ simulations. 

We generate a series of simulated sources evenly spaced 
in radial (1' spacing out to ^ 10') and azimuthal directions 
(15°) from the aim point. We simulate the Chandra High 
Resolution Mirror Assembly (HRMA) PSF using SAOSAC 
and simulate ACIS-I detection using MARX. The source 
photons are sampled from a power law model with Fph = 
1.5 and A^h = 3 x 10^" cm"^ and the background photons 
from the Markevitch 'period B' ACIS-I background dataset 
{acisi_B_i0123 _bg_evt_230301 .fits)}^ The total source pho- 
tons for a given source ranges from 5 to 1000 counts, and 
the background is integrated over 10 ksec (for < 15 count 
sources) or 20 ksec (for > 15 count sources). 

We apply wavdetect to the simulated data, using a simpli- 
fied version of XPIPE, and compare the wavdetecttd source 
positions with the true positions. Fig.|2]shows the size of such 
calculated 95% error circles (95% of the true source positions 
lie within the circle) as a function of the offset from the aim- 
point (Doffset in arcmin) and the number of net counts (cn).'^ 
Based on this result, we generate the empirical formula for 
95% error circle f err of the wavdetecttd source position. 



Pen- =0.25"+ ■ 



0.1" 



+ 0.03'' 



logi()(Cn+l) 
n 2 



1 + - 



1 



offset 



logio(Cn + 2) 



+ 0.0006 



logio(Cn+l) 

^offset 



n4 



logio(Cn + 3) 



■(5) 



The simulations are performed for ACIS-I detectors only 
and we employ the formula for both ACIS-I and ACIS-S ob- 
servations. For very small count sources (c,, <C 5), the er- 
ror estimate made by the above formula can be unreasonably 
large, so we limit the positional error as Pen- < '"95% because 
Perr ~ ^"95% for a One count source with no background. 

For matching with optical or IR sources, we use the Pen 
along with other relevant information such as boresight off- 
set, optical error circle, etc (Lavcocket al. 2005; Zhao et aQ 
I2OO5I) . For the X-ray data base, PXP generates an index table 
that cross-correlates potentially identical sources that appear 
in multiple observations of the same region. First, we correct 
the X-ray source positions using the known aspect offset for 
each observation.'* Next, we consider two sources to be po- 
tentially identical (and also set flag=157), when the distance 
(A12) between the sources is less than the quadratic sum of 
the error radius of both sources, i.e. 

Al2<(Pe',l+Pei2 + 0.7'x2)»^ (6) 

where the constant term is added to account for the accuracy 
of the absolute astrometry (0.7" for 95% error for each ob- 
servation) after the aspect offset correction'^. Note Eq. (|6j is 

''^ http://cxc.harvard.edu/chart/. 

'5 http://space.mit.edu/CXC/MARX/. 

' ^ http ://hea-www. harvard. edu/ ~maxim/ axaf/acisbg/ 

This net count (Cn) is the number reported by wavdetect and it is not the 
same as A'net which is calculated by the aperture photometi-y. 

http://cxc.harvard.edu/caVASPECT/celmon/; assume 0.6" for 90% er- 
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Net Counts (STD Dev) 

A 909 (30) 
X 91 (3) 
O 44 (3) 




4 6 
offset (arcmin) 



Fig. 2. — 95% inclusion radius (Pen) for uncertainties of wavdetectsd 
source positions by SAOSAC & MARX simulations; the data points are sim- 
ulated results for a given number of detected net counts reported by wavde- 
tect, and the solid lines are the visually chosen empirical formula. Eq. Isl . 
The dotted line represents 95% PSF radius (r<)s%, §2.3), used for the upper 
limit of Pcrr. 



used only for cross-correlating (or boresighting) X-ray to X- 
ray observations; for matching with optical or IR sources, see 
( Lavcocket al. 2005: Zhao et al. 2005). 

3. INITIAL RESULTS FOR THE ANTI-GC FIELDS 

In the following, we consider the results from the unstacked 
data. 

3.1. Source statistics 

Table |6l summarizes the initial X-ray results of the ChaM- 
Plane survey on the 15 selected (see Table 1) anti-GC obser- 
vations. The table shows the summary of source counts at 
each level of the analysis. The levels are devised to select 
valid sources and to increasingly (with level) restrict system- 
atic errors. 

The wavdetect routine detected 1028 sources in the 15 se- 
lected anti-GC observations (level 0). After careful selection 
of sources by the flag information and V&V, we found 921 
reliable sources (level 1) in these regions. To interpret the re- 
sults without concern for anomalies arising far from the aim 
point, we further limit the source list. In level 2, we use only 
data from CCDs 0, 1,2, and 3 for ACIS-I observations, and 
CCD 2, 3, 6, and 7 for ACIS-S observations. In level 3, we 
select sources within 400" of the aim point for ACIS-I ob- 
servations, and all of CCD 7 for ACIS-S observations. Note 
that level 3 uses only front-illuminated (FI) CCDs for ACIS-I 
and only one back-illuminated (BI) CCD for ACIS-S obser- 
vations. 

The last three columns in Table|6llist the number of sources 
that need special attention for further analysis: the same 
sources found in multiple observations of the same field 
(flag=157; i.e., the first 2 observations listed in Table 1), 
sources that may be piled-up (flag=147), and sources with 
their source region overlapping with neighbors' source region 
(flag=141). 

In the following, we go over the basic X-ray properties of 
the above sources. The complete source catalog as well as 
access to the X-ray and optical images and optical photometry 
data for probable IDs (see Grindlay et al 2005 for examples) 



is available on-hne. 

3.2. Source distribution 

We derive the logA^-logS' distribution to review the analy- 
sis procedure and for an initial analysis of the nature of the 
source population. In general, making an estimate of the un- 
absorbed flux for sources in the Galactic Plane fields is more 
difficult than for high-latitude sources because of the diverse 
and unknown spectral types, the higher column density and 
its strong dependence on the usually unknown distances to 
the Galactic sources. Therefore, the derivation and interpre- 
tation of the logA^-logS' distributions for ChaMPlane sources 
are not trivial. For the initial analysis, we follow a practice 
similar to the one done for ChaMP (K04), which allows us to 
verify the consistency of the analysis and to compare the re- 
sults directly with the high-latitude ChaMP results. We gen- 
erate the sky coverage and the logA^-logS' distribution in the 
Sc and He band. To simplify the analysis, we use only the 
level 3 data set. We also exclude the original target (and point 
sources in the target region if the intended target contains dif- 
fuse emission or a cluster of X-ray sources; flag=154) of each 
observation to maintain the serendipitous nature of the survey. 
For example, in the NGC 1569 field (a galaxy, Obs. ID 782), 
we consider the point sources within the 3.6' diameter^'^ circle 
from the galaxy center as a part of target, and in the 3C 123 
field (a quasar, Obs. ID 829), two sources nearby 3C 123 are 
also considered so. 

3.2.1. Sky Coverage 

To calculate the amount of sky co vered, we adopt a rela - 
tively simple but accurate procedure ("Cappelluti et al."2005'). 
First, we generate a background-only image from event files 
for a given band. Using the valid source list (level 1), we re- 
move the photons in the source regions and fill the regions 
with counts consistent with the neighboring background us- 
ing dmfilth?^ Note we keep the photons in the source region 
of the target (and sources in the target region if the target is 
extended) in the background-only image because we exclude 
these sources in counting the number of sources for logA^- 
\ogS. Second, for a given position in a CCD, using the expo- 
sure map and the background image, we estimate the count 
rate limit (r;) for detection at a given signal-to-noise ratio 
(SNR= a) by 



n = a- 



1+(0.75+A^b)1/2 max(ep) 



(7) 



where A^^ is the background count in the source region, T the 
exposure ti me, a = 3 and we u se Gehrels' approximation for 
the errors ( Gehrels et alJl98^ .^^ 

We need to convert r; to the equivalent flux in order to cal- 
culate the sky coverage as a function of flux. Note that r/ is 
scaled by the exposure map in the same way as r^gt in Eq. (|3|l. 
Therefore, one can use t he sam e rate-to-flux conversion factor 
in the previous section ( ^2.2.2> to convert r/ to the equivalent 
flux. However, there is no unique conversion from the count 
rate to un-absorbed flux because sources in the Galaxy can be 

" http://hea-www.harvard.edu/ChaMPlane/. 

Optical major diameter (D25). See lHumphrev et an<2003ri . 
^ ' http ://cxc . harvard. edu/ ciao/ahelp/dmfilth.html 

The Gehrels' approximation in Eq. is optimal for cr = 1, and for 
cr > I, Gehrels et al. 1 1986) recommands a more sophisticated formula, but 
we take the above simpler approach to be consistent with the definition of 
SNR in Eq. JSJ. 
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TABLE 6 

Summary of Source Numbers for 15 selected anti-GC observations (7 ACIS-I, 8 ACIS-S Observations.) 



Stage 



Selection Rules 



Number of Sources 



Special" 



ACIS-I ACIS-S Total multiple pile-up overlap 



level wavdetect 


Everything 


448 


580 


1028 


18x2 


3 


28 


ICVCi i OCICUL VclilQ aUUH^Cft 


nag ^ i iA, iz,A, iiui iH'U 


4^1 






18x2 


J 


20 


level 2 CCD selection 


CCD 0, 1, 2, 3 for ACIS-I Obs., 
CCD 2, 3, 6, 7 for ACIS-S Obs. 


409 


445 


854 


18x2 


2 


16 




in Be 


330 


396 










SNR > 3 


mSc 
in He 


214 
189 


316 
213 










level 3 Off-axis limit 


<400" for ACIS-I Obs., 
CCD 7 for ACIS-S Obs. 


276 


242 


518 


16x2 


2 


12 




in Be 


214 


210 










SNR > 3 


in 5c 
in He 


131 
123 


169 
110 











Notes. — See Table|5]for flag definitions. 

"Special attention required: sources found in multiple observations of the same field (flag=157), potential pile-ups (flag=147), and source 
region overlap (flag=141). No source with flag=142 (severe source region overlap) is found in these fields 

''Removes extended sources and other apparently false sources. For the case of sources near a chip boundary (flag=146), they may be valid 
sources, but we do not include them for level > 1 analysis because their position and source properties are subject to large uncertainties. 



at any distance from the Earth and the column density changes 
with the distance. As for direct comparison with the ChaMP 
results, we again assume the full column density along the 
line of the sight prov ided by the extinction models. Since 
iSchleg el et alJ ([l998|) may overestimate and 'Drimm el et al.l 
J2003h may underestimate the total extinction along the line 
of the sight, we use both models (insets in Fig.|3}- 

As for spectral models, we employ a power-law emission 
model with Fph =1.7 for the Sc band and Fph = 1 .4 for the 
He band sources as in ChaMP (K04). These assumptions 
for spectral and column density models are more relevant for 
background AGN sources than for Galactic sources. 

Using the rate-to-flux conversion factor and Eq. Q, one can 
assign the flux limit for ct = 3 at any position in the detector. 
For a given flux, the sky coverage is derived by summing up 
the detector area where the calculated flux limit is less than the 
given flux. In summation, we only include the level 3 region 
where exposure map is greater than 10% of the maximum of 
the CCD to meet the source selection criterion (flag ^ 146). 

3.2.2. logN-logS 

To be consistent with the sky coverage calculation (cr > 
3), for a given band, we count the number of sources with 
SNR> 3 for the logA^-log5 distribution except for the target 
and sources in the target region, where 

Fig Elshows the logA^-log5 distribution of the sources in the 
anti-GC fields. The solid line with the (yellow) shade is for the 
extinction model by Schlegel et al. ( 1998), and the dotted line 
with the (red) shade is for the model by Drimmel et al. ( 2003). 
The shading represents the statistical uncertainties (±1ct) of 
the distribution. The ChaMP results are shown as a thin (blue) 
solid line without any shade in the plots. Table0summarizes 
the cuiTent status of source classification (only ^ 6 - 17%) by 
optical imaging and spectrocopy for the sources in Fig.|3] 



In Fig. |3j for a given spectral and extinction model, the re- 
sults for the ACIS-I and ACIS-S observations agree within 
< 2a overall, and they exhibit better agreement (< Icr) where 
the number of sources is sufficient so that the statistical fluctu- 
ations are not too overwhelming (^o < 1 x 10"'"* erg/sec/cm^ 
in the Sc band and 5o < 3 x 10"''* erg/sec/cm^ in the He band). 

Because of the relatively poor statistics, it is premature to 
attribute differences between the ACIS-I and ACIS-S obser- 
vations to any systematic bias in the analysis. For example, in 
the case of the Sc band [Fig.|3](a) & (c)], the logA^-log5 distri- 
butions expectedly show a strong dependence on the extinc- 
tion models, which differby as much as a factor of 3 (Table^. 
For Galactic populations, the extinction depends strongly on 
the distance to the sources. In addition, the assumed power 
law emission model (Fph =1.7) may n ot pr operly describe the 
X-ray emission from these sources ( ^3.3t . For soft coronal 
sources ( ^3.3> . the assumed power law emission model likely 
underestimates the source flux and thus their contribution to 
the logA^-logS' distribution, particularly in the soft band (shift- 
ing the logA^-log5 distribution to the left in the plot). 

The unknown spectral model and the inaccurate extinction 
estimate dominate the uncertainty of the logA^-log5 distribu- 
tions in the Sc band. Nonetheless, we expect our anti-GC re- 
sults to exhibit an excess relative to the ChaMP logAf-log5 dis- 
tribution. We suspect the extinction model of iDrimmel et alJ 
(|2b03 ) underestimates the total Galactic A^h values for at least 
a subset of the fields. 

In the case of the He band [Fig. 01(b) & (d)], the logA^- 
log5 distributions for the anti-GC fields are not sensitive to 
the extinction models, and the results are consistent with 
the ChaMP logA^-logS' distribution within the statistical limit. 
This indicates that most of the strong hard sources are back- 
ground AGNs. For example, judging from the statistical un- 
certainty of the logA^-log5 distributions, we estimate that for 
S> S() = 2x10"''* erg/sec/cm^ only less than - 20% of the 
sources are Galactic in Fig.|3](b) & (d). 

It is premature to draw any conclusion from source classifi- 



10000 



1000 



muffTiH^ "I '''("" 



100 



A 




10"'= 10"" 10" 

So (erg/sec/cm^ 0.5-2 keV, PLI=1.7) 



looob' 



1000 



100 



■i,ip^fi;MMfiii|/iiiip,ii'|'ii'Bi''n|i' 



(c) ACIS-S, Sc, PLI = 1.7 



95 


19 : 


650 


26 


755 


23 


782 


34 


829 


12 


2218 


5 - 


2545 


26 : 


2796 


5 : 


Totol 


150 




10"" 10"'* 10" 

So (erg/sec/cm^ 0.5-2 keV, PLI=1.7) 



lOOOOb 



1000 



100 



10 



(b) ACIS-I, He, PLI = 1.4 



1.0 



S 0.1 




676 


14 : 


1948 


6 


2553 


18 


2787 


29 ■ 


2802 


9 


2803 


16 - 


2810 


27 : 


Totol 


119 : 



10"'= 10"" 10"'' 

. .Sp I 



10"'= lO"" 10"'^ 

So (erg/sec/cm^ 2-8 keV, PLI = 1.4) 



10000 



1000 



100 



10 



"llffMflilfllMilil"i'l' i' 



(d) ACIS-S, He, PLI = 1.4 



0.1 




10"'= 10"'" 10"'^ 



10"'= 10"'* 10"'^ 

So(erg/sec/cm^ 2-8 keV, PLI = 1.4) 



Fig. 3. — The logW-logS distributions of tlie X-ray sources in the 15 anti-GC observations (90° < 1 < 270°): (a) ACIS-I observations in the Sc band, (b) 
ACIS-I in He, (c) ACIS-S in 5c, and (d) ACIS-S in He. The solid line with the (yellow) shade is for the extinction model by Schleael et al. 1 1998), and the 
dotted line with the (red) shade is for the model bv Drimmel et alJ 120031) . The shading represents the statistical uncertainties (Ict) of the distribution. The thin 
(blue) solid line is the extra-Galactic result from ChaMP (K04). The inset shows the sky coverage: the solid line from the model by Schleael et al. 1 1998), and 
the dotted line from Drimmel et ah ;.20Q3<1 . The top right list in each plot indicates the Obs. IDs and the number of sources with SNR > 3. The targets of each 
observation are excluded. The vertical black bars abov e each plot sh ow the actual source distribution by flux assuming the model by Schleael et al. 1 1998^) (top 
row) and the flux distribution assuming the model bv iDrimmel et all t20Q3i^ (bottom rows). See the electronic ApJ for color version of the plots. 



cation (Table0or Table|8} because only ^ 10% of the sources 
in Fig. |3l (b) & (d) are identified as of this writing. How- 
ever, the above log^V-logS' results are consistent with source 
classification in Table because, among 8 sources (out of 
15 classified sources in the He band) that have S > 2x10"''* 
erg/sec/cm^ and SNR > 3 in the He band, 7 sources are 
quasars and only one is Galactic. Note that in the above 

Note the logW-logS distributions in Fig.|3](and TablelTl use level 3 data 
and exclude the target of each observation. On the other hand, Tablelsldoes 
not exclude the target of observations, but uses selection criteria different 
from those in TablelTl Tablelslcontains 10 sources (out of 25 distinct sources) 



comparison, we need a limit on S to alleviate the selec- 
tion effects arising from spectral difference between stars and 
quasars in TableQ 

Assuming the extinction model bv lSchlegel etaDiMi IS 
correct or at least does not underestimate the total Galactic A^h 
value, let us compare the results using Schlegel et alJ il998t) 
in Fig.|3] The logA^-log5 distributions in the Sc band shows a 
larger excess compared to the ChaMP result than the distribu- 
tions in the He band. Several reasons may explain these dif- 

that have S > 2 X lO"''* erg/sec/cm^ and SNR > 3 in the He band, and among 
them, 8 sources are quasars and the other two are Galactic. 
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TABLE 7 

Classification of Sources in Fig.^by Optical Imaging and 



Spectroscopy 




ACIS-I 


ACIS-S 






(a) Sc 


(h)Hc 


iOSc 


(d) 


He 


Total 


126 


119 


150 


94 




"Matching performed 


89 


76 


124 


79 




''Optical Matches 


65 (6) 


38 (5) 


71 (3) 33 


(0) 


'Spectrum Observed 


35 


23 


41 


19 




Identified as Stars 


15 


5 


24 


3 




Quasars 


8 


4 


3 


3 




Galaxy 


1 













Unclear 


11 


14 


14 


13 





Notes. — The logA'-log5 distributions in Fig. |3| exclude targets of 
each observation. 

"As of this writing, source classification with optical spectroscopies 
is performed for 4 ACIS- I Obs. and 6 ACIS-S Obs. among total 15 
Obs. See also Zhao et all i2005l) . 

''() for the number of sources with multiple matching candidates. 
'See lRoeel etalj<200^ . 



ferences. First, the X-ray emission of many Galactic sources, 
e.g. coronally active stars - that are detected in ChaMPlane 
but not in ChaMP- is relatively stronger in the soft band. Al- 
though only < 15% of the X-ray sources are classified (Ta- 
ble 0, the results suggest that, as expected, the star to quasar 
ratio in the logA^-log5 distribution of the soft band is higher 
than the same for the hard band. Second, for Galactic sources, 
the extinction and henc e the unabsorbe d flux is likely overes- 
timated by the model ofl Schlegel et alJ lTlQQSl which can sys- 
tematically shift the entire curve in Fig. |3]to the right. The 
effects of extinction are relatively larger in the soft than in the 
hard band. 

To derive the true Galactic logA^-logS' distribution, we need 
to allow for the differing source spectra and distances (and 
thus A^h), and then re-derive the flux estimate of each source 
separately. In future papers, we shall incorporate additional 
information (e.g. from optical spectra and thus reddening) for 
sources whenever available for more complete logA^-log5 dis- 
tributions. 

3.3. Source Properties by QCCD 

Fig Sshows the QCCDs of the anti-GC field sources with 
A^net > 4 in By. The deta iled description of the QCCD is 
found in 'Hon g et al.1 J2004 . In summary, the jc-axis shows 
the median of the photon energy distribution (m) and the y- 
axis shows the ratio of two quartiles of the photon energy dis- 
tribution (3 Qis/Qis)- Note that the median m is related to the 
median energy (£50%) by 



where = 0.3 and £'up= 8.0 keV in the Bx band dHong et al.l 
|S)04). The top jf-axis of the QCCD is labeled by £"50% and 
the bottom x-axis is the inverse hyperbolic tangent of m (this 
allows a convenient quantity for plotting). 



Fig. 13(a) is for ACIS-I observations, overlayed with power- 
law model grids, and Fig.|4](b) shows the same data with the 
eiTor bars and thermal Bremsstrahlung model grids. Fig. |4] 
(c) and (d) are for ACIS-S observations. The figures con- 
tain sources from level 2 data for ACIS-I, and level 3 data 
for ACIS-S observations (see Table |6jl except for two piled- 
up sources (PSR J0538H-2817 and GK Persei). Note the dif- 
ference in the grid patterns between the FI and BI CCDs, 
which is mainly due to relatively high detection efficiency 
of BI CCDs at low energies. The temporal degradation of 
the low energy efficiency also makes the grid patterns change 
(shrink). The grid patterns in Fig. |3 are averaged over the 
selected observations. 

In Fig. 13(a) and (c), one can roughly identify sources in 
a few distinct groups. In particular, the sources within the 
shaded (yellow) region appear to be too soft to be described 
by power-law models (Fph > 3-4). The rest are relatively 
well described by power law models with conventional values 
for Fph ('^ 1 - 3) but reveal a wide range of extinctions. In 
Fig.|4](b) and (d), one can see the very soft sources are better 
described by thermal models (kT < 1 keV), which indicates 
many of these sources are most likely to be stellar coronal 
emission sources (stars). The hard sources that are better de- 
scribed by the power law models are likely accreting X-ray 
sources and are predominantly background AGN but can also 
include accretion-powered compact binaries (e.g. cataclysmic 
variables; see Grindlay et al 2005). 

Fig Elshows the same QCCD plots but only with the known 
sources or optically identified ones (Rogel et al. 2005). Ta- 
ble |S] contains the complete list of the sources in Fig.|5] The 
solid circles indicate stars and the open circles indicate back- 
ground AGNs. The rest are pulsars, black hole (BH) X-ray 
binary (XB), etc. The QCCD separates relatively soft stars 
and hard quasars even with a small number of net counts (> 
40 counts in Bx). Five sources in Obs. ID 2787 are also found 
in Obs. ID 1948. Note that their quantile values (and net count 
rates) from two observations are consistent, even though three 
sources in Obs. ID 1948 have only 16.8, 25.4 and 8.4 net 
counts respectively (Table|8j. 

4. SUMMARY 

We describe the X-ray analysis procedure for X-ray pro- 
cessing of ChaMPlane survey archival data using custom de- 
veloped analysis tools, which can be useful for other Chan- 
dra analysis projects as well. The initial X-ray results from 
the analysis of the 15 selected anti-GC observations reveal a 
few distinct classes of sources. The detailed X-ray analysis 
of source distributions and source classifications for all Anti- 
center ChaMPlane fields, followed by Galactic Bulge and GC 
region fields, will be presented in subsequent papers. 
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Fig. 5. — The same QCCDs as in Fig.l^but only with the optically identified sources IRogel et all2005D or known sources (A'net > 40 in Bx)- The solid circles 
are stars and the open circles are background quasars and the rest are black hole binaries, pulsars, etc. See Tablelslfor the complete source list of the figure. 
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TABLE 8 



List of optically identified or known sources with Nnct > 40 in Bx 


Source ID 
XS 


Source Name 
CXOPS 


^50% 

(keV) 


3225/275 


• net 

in Bx 


in Bv 

' net 

(cts/ksec) 


Source Type 


Known Name*" 


CCD 0,1,2,3 (FI)forAClS-I obsen-ations 














02810B2_018 


J2358 13.4+622447 


2.14(23) 


1 11(11) 


97.7 


2.4(3) 


Quasar z=:1.29 




\J\J\J 1 \J±J^ \J L J. 


T0421 1+324724 


1.71(12) 


1.33(19) 


64.2 


3.8(5) 


Quasar z= 1.845 




0067 6B0 001 


J04220 1.6+325728 


1.39(10) 


1.16(24) 


63.5 


3.7(5) 


Quasar z=2.203 




00676B3_018 


J042123.7+324836 


1.66(14) 


1.08(17) 


57.0 


3.6(5) 


Quasar 




00676B2_001 


J042133.8+325556 


1.64(30) 


0.99(14) 


48.5 


2.7(5) 


Quasar z=2.055 




"02810B0_021 


J235813.2+623343 


2.44(42) 


0.87(20) 


42.4 


1.0(2) 


Quasar z=2.04 




00676B1_006 


J042211. 8+325604 


1.67(33) 


1.12(15) 


41.8 


2.7(5) 


Quasar z=0.65 




02787B3_002 
''01948B3_001 


J222905.2+611409 
J222905.2+611409 


2.33(06) 
2.31(13) 


0.99(02) 
0.98(06) 


1972.6 
327.4 


21.9(5) 
25.2(1.5) 


Pulsar 


PSR J2229+6114 


02787B1_005 
01948B1_002 


J222833.4+611105 
J222833.3+611105 


1.10(04) 
1.13(12) 


1.52(06) 
1.64(16) 


237.9 
46.8 


2.9(2) 
3.5(6) 


YSO K star 




02787B0_003 
'=01948B0_001 


J22300 1.4+6 11059 
J223001.3+611100 


1.35(09) 
1.62(34) 


1.56(19) 
1.41(76) 


87.6 
16.8 


1.1(1) 

1.3(4) 


early K star 




02787B3_007 


J222847.3+611214 


1.65(10) 


1.42(15) 


79.1 


1.6(2) 


G star 




'=01948B1_001 


J222847. 3+61 1214 


1.76(21) 


1.45(26) 


25.4 


1.8(4) 




02787B3_005 
'=01948B3_004 


J222853. 1+611351 
J222853.0+611351 


1.35(07) 
1.40(39) 


1.58(35) 
1.38(37) 


48.0 
8.4 


0.5(1) 
0.6(3) 


G/K? star 




CCD 7 (BI) for ACIS-S observations 














00829B7_002 


J043704.3+294013 


2.92(12) 


1 26f06) 


380.7 


14.1(8) 


Quasar z=0.218 


3C 123 


00782B7_002 


J043 125.0+645 154 


3.73(12) 


1.82(10) 


157.8 


2.7(2) 


Quasar z=0.279 




00829B7_003 


J043654.8+294018 


1.69(25) 


0.81(20) 


45.1 


1.7(3) 


Quasar 




02218B7_004 


J044909.0+450039 


1.45(10) 


1.25(30) 


47.2 


2.3(4) 


Quasar 


3C 129 


00829B7_022 


J0437 18.4+294546 


1.07(05) 


1.68(10) 


136.2 


7.4(7) 


K star 




00782B7_042 


J043003.0+645143 


1.18(05) 


1.54(12) 


133.6 


1.7(2) 


YSO M star 




00650B7_015 


J033 108.2+435750 


0.94(03) 


1.42(19) 


111.3 


2.0(2) 


dMe star 




00755B7_005 


J222545.8+653833 


0.95(06) 


1.42(16) 


75.4 


1.7(2) 


dMe star 




00782B7_010 


J043058.3+644852 


1.00(04) 


1.50(23) 


63.7 


0.8(1) 


dMe star 




00782B7_018 


J043039.0+645013 


0.78(11) 


0.86(20) 


40.2 


0.5(1) 


dMe star 




00782B7_012 


J043057.4+645048 


1.72(04) 


1.34(07) 


257.8 


3.1(2) 


Nebula 




00095B7_004 


J062244.5-002044 


1.17(09) 


1.19(12) 


99.7 


3.1(4) 


BHXB 


A0620-00 


00755B7_004 


J222552.5+653535 


1.25(10) 


0.90(14) 


83.7 


1.8(2) 


Pulsar 


B2224+65 



See Eq. The top X-axis of the QCCDs is labeled by £50%. 
Ones with the known name is the target of the observation. 

Three candidates are found for the counter part. The nearest candidate is identified as a quasar and the other two are unknown type. 
This source is not plotted in Fig.|5|for clarity. 

Note A'nct < 40 for these sources (not shown in Fig.|5}, but they are listed for comparison with the same sources found in Obs. ID 2787 with 
Wnet > 40. 



